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1. INTRODUCTION 


An empirical model was developed to infer soil moisture and surface roughness 
from radar data. The accuracy of the inversion technique is assessed by comparing soil 
moisture obtained with the inversion technique to in situ measurements. The effect of 
vegetation on the inversion is studied and a method to eliminate the areas where 
vegetation impairs the algorithm is described. 

2. EXPERIMENTAL DATA AND MODEL DESCRIPTION 


Two sources of experimental data were used in the derivation of the empirical 
model. The University of Michigan's LCX POLARSCAT is a truck-mounted network- 
analyzer-based scatterometer operating at three frequencies: 1.25, 4.75 and 9.5 GHz. The 
POLARSCAT data set (Oh et al., 1992) used in this study includes the co- and cross- 
polarized (hh, w and hv) backscatters for four surfaces. The RASAM radar developed 
by the University of Bern, Switzerland is a truck-mounted radiometer-scatterometer with 
the scatterometer system operating at six frequencies between 2.5 GHz and 11 GHz. A 
complete description of the instrument can be found in (Wegmuller, 1993). The RASAM 
data set includes measurements of the hh, w, hv and vh-polarized backscattering 
coefficients over a variety of surfaces. 


Using these two data sets, the hh-polarized and vv-polarized backscattering 
coefficients and <J^ V were empirically found to follow these two equations: 

oi = lO- 275£2 ^lO O2Sfun0 (-ytsin0) ,4 A 0 - 7 

sin 6 

o°„ = io -2 io° 46Cun 0 (kh sin 6) U X 07 

sin 6 

where 6 is the incidence angle, £ is the real part of the dielectric constant, h is the 
RMS height of the surface, k is the wave number and A is the wavelength in cm. 
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The general backscatter behavior with roughness and dielectric constant is 
similar to the trend predicted by the small perturbation model and the physical optics 
model (Ulaby et al., 1986). The backscatters decrease with increasing incidence angle 
and with decreasing roughness. Restricting the validity of the model to kh £ 2.5 and 
6 2 30 ° will insure that the ratio is always less than 1 . We note that the 

roughness of the natural surfaces that the algorithm is applied to rarely exceeds kh = 2. 5 
(corresponding for L-band to an RMS height of 10 cm). 


3. EFFECTS OF VEGETATION AND CALIBRATION 


The inversion relies on the co-polarized channels and does not use the cross- 
polarized channels. This has several advantages. The First advantage is that the co- 
polarized channels can be calibrated directly with passive targets like corner reflectors. 



The cross-polarized channel calibration is usually derived from measurements made on 
the co-polarized channels and as a result is less accurate than the calibration made on the 
co-polarized channels. The second advantage is the normally higher signal-to-noise ratio 
in the co-polarized channels than in the cross-polarized channels. Finally, the third 
advantage is that vegetation is known to have a relatively strong effect on the cross- 
polarized channel (Le Toan et al., 1992). It can therefore be expected that an inversion 
algorithm relying on the co-polarized channels will be more robust to the presence of 
vegetation than one relying on both the cross- and the co-polarized channels. 
Nevertheless, for a significant amount of vegetation, the algorithm will overestimate 
surface roughness and underestimate soil moisture. 

Active microwave sensors have been shown to be good discriminators for 
biomass (Le Toan et al., 1992). This capability can be used to select the areas with low 

vegetation cover where the inversion can be applied. TheL-band CT^/<7^ ratio image 

acquired over Chickasha, Oklahoma was compared to a SPOT-derived Normalized 
Difference Vegetation Index (NDVI) image (tucker, 1979) over the same area. Overall, 
a pixel-to-pixel comparison between the two indices provides the regression curve. This 
curve shows that the L-band parameter is positively correlated to the NDVI. We found 

that masking out the areas for which the L-band <T° V / CF^, ratio is greater than -11 dB, 

results in a reliable soil moisture estimate. The -1 1 dB O'J, /<7° ratio corresponds to 
NDVI of around 0.4. 

The study showed that a desired 4% accuracy in soil moisture estimate requires 
a .5 dB accuracy of the relative calibration and a 2 dB accuracy in the absolute 
calibration. Both requirements are met by both the AIRSAR L-band and the SIR-C 
sensor. It should therefore be possible to derive accurate soil moisture maps for the data 
provided operationally by these sensors. 

4. SAR DATA INVERSION 

The inversion technique was first applied to Washita'92 AIRSAR data. The 
experiment followed a period of heavy rain so that the conditions on June 10, 1992 were 
very wet with standing water and saturated soils fairly common. No further rain fell 
during the next nine days and we were able to follow a drying pattern. The area covered 
in the AIRSAR Workshop Slide 1 is around 8 km by 10 fan. It was imaged by the 
AIRSAR system on six different days between June 10 and June 18, 1992. An L-band 
hh image is displayed at the lower right comer. The AG002 bare field which was 
extensively studied is outlined in the image. The other six images were obtained by 
inverting the L-band data acquired on the corresponding days. The dielectric constant 
maps were translated into soil moisture maps using the Hallikainen empirical curves 
(Hallikainen et al., 1985). The black areas in the soil moisture maps indicate areas where 
the vegetation cover is too thick for the inversion to be reliable. The scatter plot is a 
comparison between the in-situ measurements and the estimated values of soil moisture 
over the AG002 field. The RMS error of the soil moisture estimate is 1.6 %. The RMS 
height maps corresponding to the six data takes are very similar and indicate no trend in 
roughness changes. The quantitative results are presented in Table 1. 

Washita^ was an aircraft and shuttle experiment that occurred between April 9 
and April 18 in which the shuttle took data on April 1 1 through 17 and the AIRSAR on 
April 10 and 11. Two of the sampled fields where data are already available are in the 
radar scene and meet the low vegetation criteria described earlier. The results from the 
inversion are included in Table 1 for both SIR-C and for AIRSAR. 

The soil moisture inversion algorithm was also applied to an AIRSAR data set 
acquired over Spain in the summer of 91 during the EFEDA campaign (Bolle et al., 
1992). Three of the sampled fields are in the radar scene and meet the low vegetation 
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criteria described earlier. The results from the inversion are included in Table 1. Over all 
the areas where SAR data and ground truth measurements are available, the RMS error in 
soil moisture is 3.5% (see Table 1 and Figure 1). 

5. CONCLUSIONS 

We presented an empirical algorithm to infer soil moisture from imaging radar 
data over bare surfaces using two co-polarized radar cross-section measurements. The 
algorithm was developed with scatterometer data, and tested with several data sets 
acquired with the AIRS AR system; we also presented the first soil moisture images 
derived from spacebome SIR-C SAR data. A comparison with in situ data shows that the 
algorithm infers soil moisture with an accuracy of better than 4 %. Best results are 
achieved when the surface roughness is such that kh < 2.5 (10 cm RMS height for L- 
band) and the incidence angle is larger than 30 degrees. We also quantified the 
calibration requirements of the algorithm and established that current operational 
multipolarization SAR systems such as AIRSAR and SIR-C routinely deliver images that 
meet or exceed these calibration requirements. We also presented evidence that the ratio 
of the cross-polarized return to the like-polarized return could be used to decide which 
areas the algorithm can be applied to and showed that this ratio was positively correlated 
with NDVI. 


Scene 

Sensor 

Data 

Take 

Date 

Field ID 

Mv 

Msd 

[%] 

Mv 

Est. 

[96] 

P| 

h 

Est. 

(cm) 

CHICKASHA 

AIRSAR 

3902 

6/10/92 

AG002 

28.7 

29.2 

1.19 

1.18 

CHICKASHA 

AIRSAR 

3664 

6/1 2/92 

AG002 

22.4 

21.2 

1.19 

1.35 

CHICKASHA 

AIRSAR 

3359 

6/13/92 

AG002 

24.1 

23.5 

1.19 

1.09 

CHICKASHA 

AIRSAR 

3875 

6/14/92 

AG002 

18.1 

19.4 

1.19 

1.19 

CHICKASHA 

AIRSAR 

3883 

6/1 6/92 

AG002 

13.6 

17 

1.19 

1.24 

CHICKASHA 

AIRSAR 

3360 

6/1 8/92 

AG002 

11.6 

12.1 

1.19 

1.45 

CHICKASHA 

AIRSAR 

4040 

6/17/92 

AG002 

17.5 

18.5 

1.19 

1.45 

CHICKASHA 

AIRSAR 

4041 

6/13/92 

AG002 

24.1 

27.4 

1.19 

1.38 

CHICKASHA 

AIRSAR 

Mosaic 

6/1 3/92 

RG148 

27.6 

30.6 



CHICKASHA 

AIRSAR 

Mosaic 

6/1 3/92 

RG131 

29.2 

34.5 



CHICKASHA 

AIRSAR 

4254 

4/12/94 

Field 12 

18.4 

24.5 


1.58 

CHICKASHA 

AIRSAR 

4254 

4/1 2/94 

Field 1 5 

24.8 

21.8 


0.86 

EFEDA 

AIRSAR 

3146 

6/13/90 

Field 2 

EQ 

6.9 

1.41 

1.34 

EFEDA 

AIRSAR 

3146 

6/13/90 

Field 4 

30.6 

22.8 

0.6 

0.9 

EFEDA 

AIRSAR 

3146 

6/1 3/90 

Field 5 

18.6 

18.2 

1.79 

0.8 

CHICKASHA 

SIR-C 

10155 

4/12/94 

Field 12 

18.4 

20.2 


1.2 

CHICKASHA 

SIR-C 

10155 

4/12/94 

Field 1 5 

24.8 

30.3 


1.24 

CHICKASHA 

SIR-C 

10158 

4/15/94 

Field 1 2 

9.9 

12.5 


2 

CHICKASHA 

SIR-C 

10158 

4/15/94 

Field 15 

12.5 

11.8 


1.5 





j 





|RMS ERROR 





3.3 


0.34 


Table 1: Comparison between in situ measurements and estimated values. Mv stands for 
volumetric soil moisture, h for the RMS height Est and Msd indicate the radar derived 
values and the ground truth measurements. 
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Figure 1: Radar derived soil moisture versus in situ measurements. 
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